Recent Res. Devel. Stat. Phys., 2 (2002) 83-106 



Multiparticle random walks 

L. Acedo and S. B. Yuste 
Departamento de Fisica, Universidad de Extremadura, E-06071 Badajoz, Spain 



ABSTRACT 

An overview is presented of recent work on some statis- 
tical problems on multiparticle random walks. We con- 
sider a Euclidean, deterministic fractal or disordered lat- 
tice and N ^ 1 independent random walkers initially 
{t = 0) placed onto the same site of the substrate. Three 
classes of problems are considered: (i) the evaluation of 
the average number {SN{t)) of distinct sites visited (ter- 
ritory explored) up to time t by the N random walkers, 
(ii) the statistical description of the first passage time 
tj^N to a given distance of the first j random walkers {or- 
der statistics of exit times), and (iii) the statistical de- 
scription of the time tj,N elapsed until the first j random 
walkers are trapped when a Euclidean lattice is randomly 
occupied by a concentration c of traps [order statistics 
of the trapping problem). Although these problems are 
very different in nature, their solutions share the same 
form of a series in ln~"(iV) In™ In(iV) (with n > 1 and 
< TO < n) for N ^ 1. These corrective terms con- 
tribute substantially to the statistical quantities even for 
relatively large values of N. 



I. INTRODUCTION 

It is now a commonplace in the history of physics to 
cite the works of the botanist Robert Brown (circa 1827) 
on the stochastic movement of pollen grains suspended 
in water as the starting point in our understanding of the 
microscopic basis of diffusion. In the "annus mirabilis" 
of 1905, Einstein gave a successful theoretical explana- 
tion of the phenomenon in terms of an atomistic theory 
which was later brilliantly confirmed by the experiments 
of Perrin and served to remove the reluctance toward the 
atomic hypothesis that scientists as eminent as Ostwald 
and Mach still entertained 0. It was in this context 
that the theory of a single random walker emerged. In 
the simplest random walk model a particle occupies a 
site of a lattice (regular, fractal or disordered) and per- 
forms a jump to a randomly selected nearest neighbour 
of that site every time step [3. Random jumps are an ef- 
fective way of simulating the net fluctuating force experi- 
mented by the Brownian particle suspended in the liquid. 
This discrete definition of the random walk has gained 
widespread acceptance since the availability of computers 
because, as are all lattice and automata models, it is es- 
pecially suitable for computer simulation. Random walks 
are also a way of describing fluctuations in the diffusion 
process that are completely smoothed out when we take 



the continuous limit represented by diffusion equations. 
Random walks have also been a very useful tool in the 
study of transport in such disordered media as fractured 
and porous rocks, silica aerogels and percolation clusters 
3], substituting more phenomenological approaches Q, 
and have also been used as topological models of poly- 
mers (the so-called self-avoiding random walks 5]). 

The single random walker statistical problems have 
been the subject of intense research since the beginning 
of past century and constitute now, in many areas, an al- 
most closed discipline extensively treated in general ref- 
erences 0. However, the generalization of these prob- 
lems to the case of iV > 1 interactin g [al or independent 
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random walkers have only really started to be consid- 
ered in detail in the last decade. It may seem strange 
to the layman in random walk theory that there are any 
reasons for such a late study of the multiparticle case, 
especially in the absence of random walker interactions, 
because in other fundamental physical models composed 
of independent particles, such as the ideal gas, the quan- 
tities of interest are obtained as simple averages over the 
mechanical properties of single particles. However, there 
are some quantities defined in a multiparticle random 
walk that can not be analyzed in terms of the single 
walker theory even when the walkers are independent! 
Examples are the number of distinct sites visited by a 
set of N independent random walkers, S^it), or the ar- 
rival time of the j-th particle of a set of N independent 
random walker at a given border, tj^N, neither of which 
yield to that simple approach. Fortunately, these prob- 
lems are tractable in the limit of large N by resorting to 
asymptotic [ll 111 111 111 Ijl O or Taube- 

rian and Abelian techniques HHHHI. The complexity 
of these problems is a consequence of the superposition 
of the trails of the random walkers in the case of the 
territory problem and the competition between random 
walkers in the case of the order statistics. Hence every 
random walker has an influence on the result which is in- 
directly correlated with the rest of the random walkers' 
influences in spite of the absence of direct interactions 
between them. 

Interest in multiparticle diffusion problems has been 
rekindled lately by the development of experimental tech- 
niques allowing the observation of events caused by single 
particles of an ensemble J2, 23J. These are powerful tech- 
niques aimed at the study of local conditions (mechanical 
response, viscoelasticity) inside such complex structures 
as fibrous polymers, the intracellular medium, etc., which 
determine the behaviour of molecular motors and the rate 
of biochemical reactions. These new research tools for 



the analysis of soft disordered media also demand a b( 
ter understanding of the statistical problems associat 
with mesoscopic samples of Brownian particles. 



II. STATEMENT OF THE PROBLEMS AND 
DEFINITIONS 

A. Substrates 

We will consider three different classes of substrates 
the following sections: (i) d-dimensional Euclidean 1; 
tices {d = 1, 2, 3), (ii) deterministic fractals (in particuli 
the two-dimensional Sierpinski gasket) and (iii) stoch; 
tic fractals (in particular, the two-dimensional incipie 
percolation aggregate embedded in the square lattice) 

Euclidean lattices are reminiscent of crystalline stri 
tures and are the most customarily studied media, 
that the solution of our problems in this case is func 
mental. Fractal lattices are good models of disorder 
substrates and have been used extensively in the ana' 
sis of diffusion and transport properties in these med 
In contrast with classical diffusion, the mean-square d 
placement of a random walker in fractal lattices is giv 
by (r^) ~ 2Di^/''"', where > 2 is the anomalous difi 
sion exponent {dw = 2 for classical diffusion) and D is t 
diffusion coefficient. This slowing down of the transpc 
has also been observed in disordered media, which cc 
sequently supports the use of fractals as models of r( 
disordered materials Deterministic fractals are cc 
structed by iteration of an unvariying rule 0, |^ , starti 
with a seed and successively applying the same itei 
tor. These kinds of fractal are only mathematical ic 
alizations but they have the obvious advantage of bei: 
suitable for the use of exact renormalization techniqi 
in the calculation of many statistical quantities \12L ^ 
Unfortunately, stochastic fractals do not allow the app 
cation of these techniques because they are the result 
some random process. However, they are also closer 
real disordered media which share their statistical-fraci 
structure HHiHl^. 

The most widely known stochastic fractal used in the 
retical and simulation studies is the percolation aggregc 
0, in 113. A percolation aggregate is a cluster of si1 
in a regular lattice that are assumed to be connected 
bonds between nearest-neighbours. This cluster is cc 
structed by the random filling with probability p of t 
nodes of a regular lattice. At a certain critical conce 
tration pc of occupied sites an infinite cluster, called t 
incipient percolation aggregate, appears. The incipie 
percolation aggregate is a fractal with fractal dimensi 
df = 91/48 if embedded in two dimensions and dj 
if constructed on a three-dimensional lattice. In our sii 
ulation we considered only the two-dimensional incipie 
percolation aggregate on the square lattice construct 
by the standard Leath method 26, 27] using the val 
p = Pc = 0.5927460 • • • corresponding to site percolati 
in the square lattice l^^. In Fig. Qwe show a portion 




FIG. 1: Territory explored (black points) by = 1000 ran- 
dom walkers diffusing on an two-dimensional percolation ag- 
gregate (gray points) after one thousand steps (t = 1000). 
The starting point is at the center of the figure. 



an incipient percolation aggregate in a two-dimensional 
square lattice with L x L sites. Multiparticle random 
walk problems with independent random walkers have 
been posed and studied on all of these media. The ones 
to be considered in this review are described below. 



B. The number of distinct sites visited. Territory 
explored 

A set of N random walkers are placed on the same site 
of a lattice at t = 0. As time runs the random walk- 
ers move independently, jumping randomly from the site 
they occupy to any of its nearest-neighbours and imprint- 
ing every site they visit. We are concerned with those 
sites visited by any of the random walkers. Successive 
visits are not relevant to register a site as visited. The 
number S'Ar(i) of the imprinted sites, that is, the number 
of distinct sites visited at time t is the territory explored 
by the diffusing random walkers. This is the magnitude 
whose statistical distribution we are interested in. 

The case N = 1 was posed at the beginning of the 
1950s by Dvoretzky and Erdos and has been thor- 
oughly studied since then 13. The problem was taken up 
again by Larralde et al. [^J^l who systematically treated 
the multiparticle [N > 1) version. In these pioneering 
works the general features of the solution in the limit 

^ 1 were unveiled, and they found three time regimes 



as follows 

{const X f^, t <^ ty, 

const X In''/^ (x) , tx < i < t'^ (1) 
N{Si{t)), t'^^t 

where x = N for d = 1, x = N/ In t for d = 2 and 
X = N/y/t for d = 3 The behaviour of the territory 

covered by a single random walker, (5i(t)), is also well 
known: (S'i(i)) ~ const xi-^/^ for d = 1, {Si{t)) ~ const x 
t/lnt for d = 2, and (S'i(i)) ~ const x i for d = 3. The 
existence of these three regimes is easily explained: 

a. Regime I. In this case there are so many parti- 
cles at every site that all the nearest neigbours of the 
already visited sites are reached at the next step. It 
is clear that the territory covered by the random walk- 
ers grows as the volume of a hypersphere of radius t, 
(Sn) ^ const X f^. The crossover time from regime I to 
regime II, tx, is simply derived if we take into account 
that regime I must break when the number of particles 
on the outer zone of the territory visited is of order 1. For 
very short times the number of particles on the outer vis- 
ited sites decreases as N/z*, where z is the coordination 
number of the lattice, and, thus, the overlapping regime 
will break approximately when iV/z* ~ 1 or, equivalently, 
when = 0{lniV). 

b. Regime II. This is the most interesting regime. 
The simple regime I has broken and the random walkers 
move diffusively so that the radius of the territory ex- 
plored grows as t^/^. Then the territory explored (S'Ar(t)) 
is given, essentially, by a volume that grows as mod- 
ified by a factor that depends on TV. The explored region 
is divided into an inner hyperspherical core and a corona 
of dendritic nature characterized by filaments created by 
the random walkers wandering in the outer regions (Fig. 
O. The overlapping of the trails of the random walkers 
make the problem non-trivial 0, In Refs. [lli fl^ . 
asymptotic techniques were used to obtain the prefactor 
and corrective terms of the main term of {SN{t)). 

c. Regime III. Particles are very far from each other 
and their trails (almost) never overlap so that (Sn) ^ 
N{Si). The territory explored by the set of N random 
walkers is the sum of the territories explored by N single 
random walkers. The crossover time from regime II to 
regime III is t'^ = O(e^) for d = 2 and t'^ = 0{N^) for 
d = 3 0, 01 . This regime never appears if the spectral 
dimension dg = 2df/dw of the substrate is d^ < 2. In 
particular, regime III is never reached in one-dimcnsional 
lattices. 

The territory problem was also studied for fractal lat- 
tices with ds < 2 by Larralde et al. Q and Havlin et 
al. %] who proposed the expression (Sn) ~ const x 
trf./2 (inTv)'^^/" for the regime II with iV > 1, where 
u = dyj/{dw — !)■ More recently, Drager and Klafter [l^ 
have also analyzed this problem using scaling arguments 
finding that (Sn) = O [f^'/^iln NY'/""] for < t, 
where v = d:^/(d^ — 1) and d:^ = dm/dmin is the chemical- 
diffusion exponent 00, US- Of course, the two predic- 




FIG. 2: A snapshot of the set of sites visited hy N = 1000 
random walkers on the two-dimensional lattice. The visited 
sites are in white, the unvisited ones are in black and the in- 
ternal gray points are the random walkers. The outer white 
circle is centered on the starting point of the random walk- 
ers and its radius is the maximum distance from that point 
reached by any walker at the time the snapshot was taken. 
The internal black circle is concentric with the former but 
its radius is the distance between the origin and the nearest 
unvisited site. 

tions agree for those media, such as Sierpinski gaskets, 
for which d,nin = 1 but disagree for those with dmin ^ 1, 
such as the two-dimensional and three-dimensional in- 
cipient percolation aggregates for which dmin — 1.15 and 
dmin — 1.33, respectively. It is notable that both Havlin 
et al. and Drager and Klafter supported their conjectures 
by comparison with simulation results obtained for two- 
and three-dimensional percolation aggregates. We have 
shown that their collapsing plots are not conclusive be- 
cause in these plots the influence of the large logarithmic 
corrective terms is not considered llTl. 



C. Order statistics of exit times 

Order statistics is a relatively young field of mathe- 
matical statistics j^^. Its objective is the ordering of 
sequences of a set of random variables. Ours is also 
a problem of order statistics for the sequence of first- 
passage times (exit times) of a set of N independent ran- 
dom walkers, all starting from the same site at the same 
time, to a "spherical" boundary or radius z in Euclidean 
and fractal media. In particular, we are interested in the 
calculation of the moments {tj^N{z)) of the jth passage 
time tj,f^(z). This is the time taken by the jth random 
walker of a set of N to first reach a given distance z. 
Some early results concerning the order statistics of a set 
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of random walkers on Euclidean lattices were obtained 
by Lindenberg et al. 30|J and Weiss et al. Asymp- 
totic results {N ^ 1) for finitely ramified fractals were 
obtained much later by one of us |l2| using renormaliza- 
tion techniques developed by van den Broeck ,24]. This 
kind of fractal has the property that, by cutting a finite 
number of bonds, a portion of the lattice with certain 
number of generations becomes isolated from the rest of 
the lattice. Curiously, rigorous asymptotic expansions 
for Euclidean lattices of arbitrary dimension were only 
given later by Yuste et al ^3 • 



D. Order statistics of the trapping problem. 

The "trapping" problem has for decades been one of 
most widely studied the areas of random walk theory 
[2', . In its simpler version we have a Euclidean lattice 
where a population of traps occupies the sites with prob- 
ability c. A single random walker starts moving from an 
empty site until it is absorbed by one of the traps. The 
statistical quantity of interest is the survival probability 
<i>i(i) that the random walker is not trapped by time t. 
This problem has its origin in Smoluchowski's theory of 
coagulation of colloidal particles ^ ^ and has been 
applied to many systems in physics and chemistry such 
as trapping of mobile defects in crystals with point sinks 
132, 33, 34], the kinetics of luminescent organic materials 
|33j] , anchoring of polymers by chemically active sites |35| 
and atomic diffusion in glasslike materials [l^, among 
others. The generalization of the trapping problem to TV 
independent random walkers was only considered very 
recently by Krapivsky and Redner [19l ] who studied a 
predator-prey problem in which a static prey or "lamb" 
is captured by one of a set of N diffusing predators or 
"pride or lions" in one dimension. In their problem the 
N predators are placed initially at a given distance from 
the prey. Later, the case was studied of a stochastic dis- 
tribution of prey on the half-line 20] . Also very recently 
PH an approach has been made to the evaluation of the 
rath moment (iJV) of the time tj^N elapsed until the 
first j random walkers of a set of N are trapped in a 
d-dimensional Euclidean lattice populated by traps at a 
concentration c. In Sec. I VII we will discuss the solution 
to this problem. 



III. SURVIVAL PROBABILITY: ABSORBING 
TRAPS AND BOUNDARIES 

The calculation of the average territory covered by 
iV ^ 1 random walkers and the order statistics of exit 
and trapping times requires the previous knowledge of 
the swrmwa/ functions, or their complementaries, the mor- 
tality functions for short times. Suppose we have a cir- 
cular boundary of radius r in a lattice and a single ran- 
dom walker starting at the center of this circle at t = 0. 
The probability that this diffusing particle has reached 



the distance z during the time interval (0, t) is called 
the mortality function. We will denote it as hsiz^t) 
where the subscript B indicates that this quantity refers 
to the boundary and z denotes any distance defined on 
the substrate. In stochastic fractal lattices it is com- 
mon to define the chemical distance £ as the shortest 
path measured along the lattice bonds 3]. Between this 
chemical distance and the Euclidean distance there ex- 
ists a scaling relation, £ const x r''""", with dmm ~ 1 
for Euclidean and deterministic fractals (Sierpinski gas- 
ket, Given-Mandelbrot curve, etc. . . ) and d^nin > 1 for 
stochastic fractals (for example, dmin = 1.14(2) for the 
two-dimensional incipient percolation aggregate). Thus, 
in the case of stochastic fractal lattices there are at least 
two independent ways of defining a "spherical" bound- 
ary as the set of sites with either constant £ (chemical 
boundary) or constant r (Euclidean boundary) from the 
origin. The anomalous diffusion coefhcient = d^/dmin 
appears in Einsteins relation, (z^) ~ 2Dt^/'^^^ for the 
average square distance traveled by the random walker 
by time t, so that = d„ for z = r and = d^ 
for z — £. The complementary of hB{r,t) is the sur- 
vival probability in the case of an absorbing boundary: 
T Bir^t) = 1 — /iB(r, <). Similarly, we will define the mor- 
tality function of a single random walker starting from 
the origin site of a lattice with a trap located on site r 
as the probability that site r has been visited by that 
random walker in the time interval (0, i). We will denote 
it as hp(y,t) in order to distinguish it from hB{r,t). 

These functions have been studied through the appli- 
cation of various techniques: renormalization in the case 
of finitely ramified fractals 0, |3i 113' solution of 
the diffusion equation for Euclidean media 0, and 
computer simulations [T^ l39l ] . The resulting expression 
in the short-time (large-distance) limit is given by 

hB{z, t) ^ Ar^"e-^«" {1 + /lie-" + •••}, (2) 

with e = zjlz^YI"^ = z/(\/2i5ti/<) > 1 and v = 
d^/{dlj — 1). The parameters d^. A, fi, c and hi are 
listed in Table for several lattices. 

It is known that the mortality function of a single ran- 
dom walker starting at a distance z from a lattice trap site 
is given by a an equation formally identical to |(2Jl in the 
case of Euclidean lattices [18, 38] for ^ ^ 1- This form 
also encompasses the cases of deterministic and stochas- 
tic fractal lattices, as simulation results have shown ^3 ■ 
Hence, in general, we will write 

hp{z,t)^ Ar^^'e-'^^ [1 + hir + , (3) 

for e 3> 1, and where A, /t, c and hi is a set of parame- 
ters characteristic of the lattice. In Table Hll the values of 
these parameters are listed for several lattices. In the ab- 
sence of exact theoretical expressions for hp{z,t), these 
parameters are obtained from comparison with simula- 
tion results. Figure 131 shows an example. In this figure 
the theoretical short-time behaviour of hp{z, t) given by 
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TABLE I: Parameters appearing in the asymptotic expression of the mortality function hB{z,t) of a random walker starting 
at the origin with a boundary of trapping sites at a given distance, Eq. ((2J for four substrates: the symbol dSC refers to the 
d-dimensional simple cubic lattice, Sd to the d-dimensional Sierpinsky lattice, GM to the Given-Mandelbrot curve (d^ — d^ 
for these cases) and A2E (A2C) to the two-dimensional incipient percolation aggregate with Euclidean (chemical) boundary. 
These parameters are analytical for Euclidean lattices [iM l3q|. numerical for the deterministic fractal lattices [l^. |40| and the 
output of a numerical fit to simulation results for the percolation aggregate I39II . In this case a 2000 aggregates average and a 
time interval [0, 1000] were used. 



Case 




2D 


A 


M 


c 


/ii 


dSC 


2 


1 


2{dl2Y'''-^ 
V{d/2) 


1 - d/2 


d/2 


d- 3 
2d 


S2 


In 5/ In 2 


1.05 


2.46 


1/2 


0.981 


-0.56 


S3 


In 6/ In 2 




3.36 


1/2 


1.31 


-0.46 


GM 


In 22/ In 3 




2.5 


1/2 


1.10 


-0.6 


A2E 


2.8 


1.4 


1.6 


-1.8 


2.1 




A2C 


2.4 


1.2 


1.1 


-0.4 


1.1 






|ln h\ 

In^ 

FIG. 4: Plot of the histogram J\f{\nhp) versus |ln/ip| for 
FIG. 3; Plot of ln[— \nhp{l, t)] versus In ^ averaged over 2000 the two-dimensional incipient percolation aggregate for fixed 
two-dimensional incipient percolation clusters. The trap was r and t (sohd hue) and for fixed I and t (dashed hue). The 
always placed at a site a distance ^ = 80 from the origin. The values are r = 30, ^ = 80, and t = 1000. 
line represents the function of Eq. ^ with the parameters 
listed in Table lUl 



Eq. © is compared to simulation results for the two- 
dimensional incipient percolation aggregate. One sees 
that the values of d^, A, fi, c listed in Table UTI lead to 
good agreement between Eq. (0) and the simulation re- 
sults. 

Finally, it is interesting to note that, for disordered 
media, the territory explored and the order statistics of 
the exit times are better described in terms of the chem- 
ical distance 0, 0, [s^ than in terms of the Euclidean 
distance. The reason for the advantage of the chemical 
distance description is to be found in the broadness of the 
distribution of the mortality function at fixed time t and 
Euclidean distance r in comparison with the correspond- 



ing distribution for fixed chemical distance £. Figure 01 
shows the histogram for the values of hp{z,t) for the 
values r = 30 and ^ = 80 at i = 1000 plotted from the 
results for 2000 percolation clusters. The same Euclidean 
distance in a disordered lattice corresponds to very dif- 
ferent chemical paths depending on the holes that may 
block this path between the origin and the destination 
sites. Hence, the mortality function should exhibit large 
deviations from one realization of the fractal lattice to 
another. 
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TABLE II: The same as Table^Jbut for the mortahty function of a single random walker starting at a distance z from a trapping 
site, hp{z,t). Results for the two-dimensional incipient percolation aggregate with a trap at a given Euclidean distance from 
the origin are not shown because the broadness of the distribution does not allow any reasonable fit for the limited number of 
aggregates used in our simulations (See Fig.HJ. The parameter p is [2(2D7r)^/3]^/^p(0, 1)], where p(0, 1) ~ 1.516386 Ql- The 
chemical dimension de and the volume Vb of a hypersphere of chemical radius i = 1, are also listed. 



Case 




de 


Vo 


A 


A 


c 


hi 


ISC 


2 


1 


2 




1/2 


1/2 


-1 


2SC 


2 


2 


TT 


l/lnt 


1 


1 


-1 


3SC 


2 


3 


37r/2 


1/pVi 


1 


3/2 


-1/3 


S2 


In 5/ In 2 


In 3/ In 2 


3 


0.61 


1/2 


0.98 


-0.56 


A2C 


2.4 


1.65 


1.1 


1.0 


0.8 


1.05 





IV. TERRITORY EXPLORED 

A. Asymptotic expressions 

In this section we give the main and two first correc- 
tive terms of the asymptotic expression (a power series 
in In^^ A^) of the territory explored (S'Ar(t)) by iV 1 
random walkers, thus going beyond the leading terms 
(without prefactor) discussed in Sec. ^ In the previous 
section we have defined the survival probability Tp{r,t) 
as the probability that a site r has not been visited by 
a single random walker by time t. The average number 
of distinct sites visited by N independent random walk- 
ers, (S'ivCi)), is then simply related to rp{r,t) as follows 







N 



5:{i-[rp(r,tr}) 



(4) 



where the sum runs over all the sites of the lattice. 
One must notice that there are two averages implicit 
in Eq. Q): (a) an average over all exploration experi- 
ments performed on the same lattice, represented by the 
sum X) |l ^ [rp(r,<)]^|; and (b) a second average (. . .) 
over all possible stochastic lattices compatible with the 
generation rules. Of course, in the case of determinis- 
tic lattices (Euclidean, Sierpinski gaskets, etc. . . ) only 
the first average is necessary. As the histogram of the 
survival probability Tp(z,t) is very narrow in stochastic 
fractal lattices when the trap is placed at fixed chemical 
distance z — i (see Fig. 0J , it is not difficult to see that 
Eq. can be approximated by 



N 



Y,{l-[rp{i^,t)f}{n{m)) , (5) 



where {n{m)) is the average number of sites separated 
from the origin by a chemical distance in the range 
[£„, = mAe,£rn + M], m = 0, 1,2,--- Evidently, Eq. 
(O is exact when applied to Euclidean lattices and de- 
terministic fractals and, for the sake of generality, we will 
take it as the starting point of the subsequent derivations 
for all classes of substrates. 



First, we replace Eq. jSJ by its continuum approxima- 
tion 



{Sn)^ / {l-[Tp{i,t)f}d,V„£'''-h 



(6) 



where de is the chemical dimension of the fractal lat- 
tice and dV{£) = Vq dei'^''~^d£ is the average number of 
fractal sites placed at a chemical distance between i and 
£ + d£ (values of Vo for several lattices are given in table 
Inj) . The asymptotic evaluation technique for (Sn) is in- 
spired in the behaviour of 1 — [rp(f , t)] for a fixed time 
t. This function is plotted in Fig. \E\ for several values 
of N in the case of the two-dimensional incipient perco- 
lation aggregate. We observe that it approaches a unit 
step function Q{£ — £x) when N ^ oo, £x being a value 
that depends on N. For large N, [Tp{£,t)]^ is only non- 
negligible when Tp{£, t) is very close to 1. This occurs in 
the limit ^ = £/L{t) > 1 {L{t) - V^Dt'^^'^'^ is the root- 
mean square chemical distance traveled by a single ran- 
dom walker by time t). On the other hand, l—[Tp{£, t)]^ 
approaches rapidly the value 1 as ^ decreases (large times 
and short-distance limit). Therefore, it is clear that 
l — [Tp{£,t)]^ Ri Q{£ — £x) is a reasonable approximation 
which improves as N increases. The threshold chemical 
distance £x can be reasonably defined as the distance at 
which 1 — [rp{£,t)]^ takes the value 1/2, and from Eq. 
® we find that 1/2 w nA^;;^" exp{-c^l) and the fol- 
lowing approximation c^^ wlnTV— /iwln^x+ln2Afor^x 
ensues. We then get c£,l « In iV - /i In In iV -|- In Ac^" +\ii2 
or 

£, « i2Dr/v/< (^^y^ 

1 -/i In In N + In Ac^" + In 2 ^ 




In AT 



(7) 



The integration in Eq. lO, evaluated using the Heaviside 
step function approximation for the quantity inside the 
bracets, trivially yields {3^(1)) « Vo£'^, and after insert- 
ing the expression for in Eq. ([TJ) we finally arrive at 
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FIG. 5; Function 1 - [ri(£)]^ versus ^ = e./L{t) for (from 
left to right) = 1,10,100,1000 and 10000 where Vt(t) = 
1 - C"''"exp(-C''), V = 1.7, /i = (solid line) and = 0.8 
(dashed line). We have not plotted the unphysical values that 
appear in the case with ^ = 0.8 when ^ goes to zero. Notice 
the large influence of the subdominant power term on 
the value of [rt(^)]^ which will be reflected in the value of 
S'jv(t). The triangles mark the value of {InN/cY^^ ~ for 
iV = 10, 100, 1000, and 10000. 



FIG. 6: S = [SN{t)/VQf^'^ /{4,Dt\nN) versus l/\uN for, 
from top to bottom, dimension 1, 2 and 3 and t — 200 (inside 
time regime II). We have used A'' = 2™ with m — 3, - ■ ■ ,14 for 
d = 2, 3, and m = 3, • • • ,30 for d — 1. The numerical results 
are plotted as circles and the broken [solid] lines correspond 
to the theoretical predictions for Sjv(f) to first [second] order 
as given by Eq. I^J ■ Notice that the approximation of order 
would be a horizontal line (not shown here) passing through 
1/d. The crosses correspond to the prediction of Sastry and 
Agmon [Eq. (22) of Ref. 41| with q = 1]. The dotted lines 
correspond to the result of Larralde et al. 0) after correcting 
the amplitude of the main term (see fisf '). 



the following approximation for {S 



InN 



dgjv 



I ln2 + lnAc^ -/tlnlniV I 



The elementary approach leading to Eq. © is not exact, 
but the dominant behaviour and the form of the first cor- 
rective term found (except for the In 2 in the numerator) 
coincide with the prediction of a systematic and rigor- 
ous improvement of the approach discussed above. This 
rigorous analysis for (S'at) yields [T^IITtI] 



^ ™ (IniV)" 

n=l m=0 ^ ' 



(9) 



with 



InN 



dgjv 



(10) 



and 



,(2) 



,(2) 



.(2) 



TT 

12 



~2 



(/3 - 1)^2 . 



(11) 
(12) 

— {chi — pLUj) (13) 

(14) 
(15) 



Here oj = j + \nAc^, 7 ~ 0.577215 is the Euler constant, 
and P — d(/v = di{dl^ — l)/d^. Inserting Hll|) and the 
definition of uo in Eq. © we get an expression for {SN{t)) 
that almost coincides with that obtained before in Eq. © 
by a much simpler analysis. The only difference between 
the two expressions is that, in the first-order correction, 
the term In 2 plays in Eq. ^ the role of the Euler con- 
stant 7 in Eq. 

In Fig. © the first and the second order approxima- 
tions for (S'Ar(t)) are compared with Monte Carlo simu- 
lation results for iV = 2^, 2-*, . . . , 21^^ at t = 200 in the 
simple cubic Euclidean lattice with d = 1, 2 and 3. An 
excellent agreement was found for N > 1000. The un- 
certainty in the values of the parameters A, p., c and hi 
appearing in the general expression of the point mortal- 
ity function for fractal media in Eq. (^J does not allow a 
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FIG. 7; Plot of S'jv(t)/t'*«''''"' versus In TV in the 
two-dimensional incipient percolation aggregate for A'' — 
2" ,2^ , . . . ,2^^ . The circles [triangles] are the simulation re- 
sults for t — 1000 [t = 500] averaging over 2000 aggregate 
realizations. The dashed line is the zeroth-order theoretical 
prediction with c — 1.05 and v — 1.7 and the solid line is the 
first-order approximation with c — 1.05, v = 1.7, /i — 0.8 and 
A=l. 



clear comparison in this case. However, using the tenta- 
tive values listed in Table ITU for the two-dimensional per- 
colation aggregate, good agreement has been found with 
Monte Carlo simulation results 17] as shown in Fig. [7| 

A remarkable fact associated with the asymptotic se- 
ries expression Q for (Sn) is the large value of the 
corrective terms even for very large number of parti- 
cles TV. We have, for example. In In TV/ In TV ~ 0.075 
for one mol of random walkers TV ~ 10^"^. Hence, these 
corrective terms must be included in order to make a 
sensible estimate of the territory expected to be cov- 
ered by any number of random walkers in simulations 
or experiments. As discussed in Sec. m there are three 
time regimes in the territory problem and Eq. ^ is 
derived only for regime H. Once the trails of the ran- 
dom walkers no longer overlap (regime HI), the series 
in Eq. @ fails to converge, and the convergence con- 
dition s^'^ (In In TV)™ / (In TV)" < 1 determines the span 
of the regime II. We know that in the two-dimensional 
square lattice the parameter A = 1/ In t (see Table 1^) 
and, consequently, the first corrective term in the series 
is comparable with the main term if ] In A] ^ In TV. This 
condition implies a crossover time Tx ~ from regime 
II to regime HI. Similarly, the crossover time Tx TV^ is 
found for the three-dimensional simple cubic lattice. For 
the one-dimensional lattice and any fractal lattice with 
spectral dimension < 1 the parameter A is time in- 
dependent and regime III is never reached as we have 
already discussed in Sect. ^ Thus a unified and purely 



analytical criteria for the time scale corresponding to the 
transition between the two regimes II and HI emerges 
naturally from the series in Eq. 0. 



B. Geometry of the territory explored 

The asymptotic series in Eq. can be used to un- 
veil the geometric properties of the set of visited sites on 
Euclidean lattices . In Fig. [21 one discerns a com- 
pact circular core and a dendritic ring composed of the 
trails of those random walkers that have travelled fur- 
ther than the rest. We define i?o(TV, i) as the average 
distance between the origin and the nearest unvisited 
site at time t (the radius of the inner compact core) and 
i?+(TV, <) as the average maximum distance reached by 
any of the TV random walkers by time t. Assuming com- 
pact exploration in the sense of de Gennes ^ 
sonable estimate of RQ{N,t) is given by [(5Ar(t))/Vo]^^'^. 
Qualitative arguments also yield an approximation for 
R+{N,t) on the d-dimensional Euclidean lattice [l^- In 
one dimension, the territory covered by TV independent 
random walkers is obviously a segment stretching from 
the maximum span on the left side of the origin to the 
maximum span on the right side and, consequently, we 
have R+{N,t) = {SN{t))\d=i/2 w [4J9HnTV]^/^ The 
d-dimensional walk of a random walker over the sim- 
ple square or cubic lattice is now decomposed into d 

orthogonal one-dimensional random walks and we get 

1/2 

R+{N,t) « [4i:>(t/d)lnTV]"'" because, on average, the 
random walker travels along each direction only the d-th 
part of the time. Taking into account that on Euclidean 
lattices di = d, df^ = v = 2 and c = d/2 as shown in 
Table IhI we can write the main term of Eq. ^ as Sn ~ 
Vq [i?+(TV, t)] . Therefore, the thickness of the dendritic 

Sn/Vo 



layer R+{N,t) - Ro{N,t) « [{SN)/Vof'' 
is given by the corrective terms in Eq. as follows: 
i?+(TV,i) - i?o(TV,t) w -(lni(i)/lnTV)i?+(TV,i). This 
means that the dendritic layer's thickness is a fraction of 
the compact core radius that grows with time as Inlnt 
for d = 2 and \nVi for d = 3 (see Table |H| for A). In 
two dimensions this ratio grows so slowly that the set of 
distinct sites visited scales as Vi, and is almost statis- 
tically self-similar in regime II as Fig. |H1 shows. For the 
crossover times Tx = 0{e^) {d = 2) and Tx = 0{N^) 
(d = 3), the dendritic ring outruns the inner compact 
core and we enter into regime HI. 



EXIT TIMES 



Now we turn to the order statistics problem of eval- 
uating the quantities (tj.iv), j — 1,2, . . ., i.e., the aver- 
age time (exit time) taken by the j-th random walker 
of a set of TV to arrive at a given border containing the 
origin site, and its moments. Specific techniques have 



FIG. 8: Four successive scaled snapshots of the set of sites visited by A'" = 700 random walkers on the two-dimensional lattice 
for times (from left to right) t = 2000, t — 4000, t = 6000 and t — 8000. The second snapshot has been shrunk by the factor 
l/\/2, the third by the factor l/\/3 and the last by the factor 1/2. 



been developed for the calculation of these quantities 
[H m El iHl m. They wiU be discussed below, but 
first it is interesting to explore the connection between 
the first passage time and the territory problems. We 
have found above that for Euclidean media and also for 
media with spectral dimension dg — 2df/dw < 2, the 
random walkers perform an almost compact exploration 
of the lattice in the sense of de Gennes . This means 
the following: If the maximum distance £ from the origin 
site is reached by any of the TV random walkers at time 
t, all the sites at a distance smaller than £ have already 
been visited. The territory explored by the N random 
walkers is roughly a hypersphere of radius £, and we can 
write: 



(16) 



Taking into account Eqs. ©, (|10|) and Hll|l one easily 
finds that Eq. (|16|l implies 




a+1 



VlniV/ 

/i In In A'^ — 7 — In Ao 



In TV 



(17) 



where a = df^ — 1 and Ao = Ac^. It is remarkable that 
this equation, derived following simple arguments, is ex- 
act to first order, i.e., all terms in Eq. 1171) are exact. A 
rigorous approach leads to the result |18| |39j : 



Nl 



' £ ' 


m(Q + l) 


C 


ma ^ 


.V2D. 




[ln(AoiV)J 





ma {fi In In AoiV — 7) 
InAo^ 



1- 



ma 



(1 + ma) ( — 



21n^(AoiV) 
7- ) + 2^7 - 2hic 
2^ (^ + (1 + ma)-f ) In In XqN 
(1 + ma )^2 In^ In XqN 
In^ \nXnN 



-O 



In' XoN 



(18) 



If we write the series in Eq. (|18|l in terms of In iV instead 
of In Ao we find that it coincides, to first order, with 
the approximation in Eq. H17|l . The moments of {t™pf) 
for j > 1 have also been calculated 0, 0, and the 
result is 



-| ma+l / J 

/2Dln(Ao7V) n=i 



— (19) 



where j — 2, 3, 

A„ = 1 + 



and 



1 



In AoiV 



(-ir 



5„(2) 
(n-l)! 



/ilnln AqA'^ 



1 



■ 7 



W- InAoiV 
In^ AoiV 



(20) 



The quantities S'„(i) are the Stirling numbers of the 
first kind 'S]. The expressions (|18|l and l|19() are also 
valid for Euclidean boundaries with minor and obvious 
changes (d^ ^ d^, € ^ r, . . . ). In Euclidean lattices 
and deterministic fractals there is no difference between 
using chemical or Euclidean distances. In Fig. ex- 
act scaled results for the N dependence of (i^^v)^^™ are 
plotted against the predictions of Eq. H18() for the one- 
dimensional lattice, the two-dimensional square lattice 
and the three-dimensional simple cubic lattice. Agree- 
ment is also very good for disordered media such as 
the two-dimensional incipient percolation aggregate, as 
shown in Fig. (^UJ. 



VI. TRAPPING TIMES 

A. Lifetime of the first trapped particle 

The "trapping" problem is a fundamental problem of 
random walk theory with a long tradition 0, [s]. In 
this problem a lattice is randomly filled with static traps 
placed at certain sites, and a particle performs a random 
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FIG. 9: Scaled moments T{m) = (2d/r2) ln(AoA''){t™Ar)^''™ 
for the first arrival of the first of A'^ particles at a prescribed 
boundary as a function of 1/lnA'^. First panel: first moment, 
m = 1 (i. e., mean first passage time); second panel: second 
moment, m — 2; third panel: third moment, m = 3. Numeri- 
cal results are indicated by circles (d = 1), triangles (d = 2), 
and squares (d = 3). Asymptotic results (cf. Eq. (1181 ') to ze- 
roth order, first order, and second order are shown by dotted 
curves, dashed curves, and solid curves respectively. 



walk on the lattice until it arrives at a site occupied by a 
trap, when it is absorbed. Most works have analyzed the 
survival probability of the random walker, defined 
as the probability of the random walker not being ab- 
sorbed in the time interval [0,t] ||45|. This problem have 
also been generalized to iV > 1 random walkers, and the 
survival probability of the set of A'^ random walkers has 
been defined analogously In this spirit, we define 




FIG. 10: Plot of r 



versus In A'' (A*' 



2 , 2 , . . . , 2 ) for the two-dimensional incipient percolation 
aggregate where 5 = l/(d^ — 1). The circles are simulation 
results for ^ = 50 and the lines are the zeroth-order (broken 
line) and first-order (solid line) asymptotic approximations 
with c = 0.9, V = 1.714 and p. = -0.4. 



the survival probability <f> n (t) that no particle of the ini- 
tial set of A^ diffusing particles has been trapped by time 
t. This survival function is given by 



$Ar(t) = ((l-c)^~(*)) , 



(21) 



where c is the concentration of traps and the angle 
brackets denote an average over all realizations of the 
A'' particle random walk on the lattice. Usually, $Ar(t) 
has been estimated using the Rosenstock approximation 
0, IE HOj in which the average of the exponential 

is identified with the exponential of the average: 



(22) 



with A = — In(l — c). Higher order terms are calculated 
in the extended Rosenstock approximation (or truncated 
cumulant expansion) first proposed by Zumofen and Blu- 
men [45j |. The first-order Rosenstock approximation is 

m 



exp 



-A (Sn) 1 + 



A Var(5jv) 
2 (Sn) 



(23) 



This means that the error made by using the zeroth- 
order Rosenstock approximation is C'(A^Var(S'jv)). Con- 
sequently, this approximation performs well if the con- 
dition A'^Var(S'Ar) ^ 1 is satisfied. In Ref. 21j, it was 
shown that Var(S'Ar) = O {f^ilnNy-^) on the simple 
Euclidean c?-dimensional lattice, and we have that the 
zeroth-order Rosenstock approximation works well when 
A^i'^(ln A^)'*"2 < 1. Hence, the approximation be- 
comes poorer as A^ increases for d = 3, and also for long 
times when, eventually, the Donsker-Varadhan regime 
settles in Hg, $Ar(i) - exp(-t''/(''+2)). If we now de- 
fine hi,M{t) = —d^N(t)/dt as the probability that the 
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FIG. 11: The function 10^/{ti,jv) versus IniV for the two- 
dimensional lattice with c = 4 x 10"" and = 2^, 2'', . . . , 2^*^. 
The simulation results are averaged over 10^ configurations 
(circles). The dotted line represents the asymptotic approxi- 
mation of order 0. We also plot the results obtained by numer- 
ically integrating Eq. (1241 when the first-order (dashed line) 
and second-order (solid line) asymptotic approximations for 
(Siv(t)) are used. 



first absorbed particle of the initial set of N disappears 
during the time interval (t,t + dt], it is clear that, using 
Eq. 122|l . the average time (lifetime) at which the first 
particle is trapped is given by 



t"'hi^N{t)dt 



'0 

m 



-cxp[-X{SN)]dt . (24) 



Inserting the main asymptotic term of (Sn) given by Eq. 
(|10|l into Eq. for the number of distinct sites visited 
into Eq. H24|l a zeroth-order approximation for the mth 
moment of the first trapping time, ti^jv was found in Ref. 
[21.] as follows: 



r(l + 2Tn/d) 



{Wo 



.2m/d 



(4D IniV)' 



(25) 



No further terms in the expansion of (Sn) were con- 
sidered because, in the two- and three-dimensional Eu- 
clidean lattice case, they depend on time, and the analyti- 
cal integration in Eq. (|24l) cannot be explicitly performed. 
In Fig. 111! numerical simulation results for (tijv) on the 
two-dimensional square lattice with c = 4 x lO"** are 
compared with the predictions of Eq. (|25|l and those of 
Eq. H24|) obtained by numerical integration using the first 
and second order corrective terms of (Sn) in Eq. I^. 



B. Lifetime of the jth trapped particle 

Until now we have discussed the statistics of the ab- 
sorption of the first particle, but the problem of the order 
statistics have also been solved for 1 < j <C iV . This 
has been possible for independent random walkers. If we 
define ^(t) as the survival probability of a single ran- 
dom walker in a given trapping configuration, then the 
distribution for the absorption of j particles is given by 



(1 - *)^ 

(^^lynf j\^N-j+m_ (26) 



E 

m— 



An average over all different trap configurations yields 
^jy{t) = (*^(t)), the survival probabihty of the N ran- 
dom walkers and $j,Ar(i) = {'^j^N{t)), the probability 
that exactly j random walkers from the set of N have 
been trapped in the time interval [0,t]. From Eq. 126(1 
we get 

'^j.Nit) - (-l)^"(^^)A^$Ar(t) (27) 



(28) 



using the backward difference formula for the jth deriva- 
tive and the approximation ~ d^ /dN^ when j ^ N. 
In Fig. ((12|l the jth survival probability ^j.N{t) ob- 
tained from Eq. H27|l using the Rosenstock approxima- 
tion for ^N{t) is compared with simulation results for 
the two-dimensional lattice with a concentration of traps 
c = 4x 10"''. 

The probability that the jth absorbed particle of the 
initial set of TV disappears during the time interval {t,t + 
dt], hj^^lt), satisfies the recurrence relation: 



(29) 



with hQ pf{t). The mth moment of the time at which the 
jth particle is trapped is given by 



''J.N/ 



t"'hj^Nit)dt 



(30) 



and, taking into account Eqs. (|29|l and (|27|l . we finally 
arrive at a recurrence relation for these moments 

{tT+i,N) = (t ■» + (-1)^' ^ ^ A,(tr:^) . (31) 

For large N and small j, this relation yields 

d"T(l-h2m/d) (IniV)-!-™ 



(32) 
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FIG. 12: The jth survival probability $j,]v versus time t for 
(from top to bottom) j = 0,1,2,3,4, with TV = 1000 and 
c = 4 X 10"'' for the two-dimensional lattice. The lines repre- 
sent i.e., the zeroth-order Rosenstock approximation 
with {5'jv(t)) given by the second-order asymptotic approxi- 
mation. The circles are simulation results averaged over 10^ 
configurations. Inset: $o,Jv(i). 



where the difference operator has been approxi- 
mated by /dN^ and (t™^) was given by Eq. (|^ . 
It is remarkable that the main asymptotic term of 
the ratio (^j.N / (t j.jv) between the variance (Tj,Af — 
\2 and (t, 



(t2 



- t 



depends only on the dimen- 



sion of the lattice as is easily shown from Eq. H32() : 



a,.N [r(l + 4/rf)-r2(l + 2/d)] 



1/2 



r(i + 2/d) 



(33) 



In Fig. El this ratio is plotted for d = 1, 2 and 3 for 
several values of j and N and compared with simulation 
results. Good agreement is found for N > 1000. 
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FIG. 13: The ratio cr^, jv/{tj,jv) , j = 1 (circles) j = 2 
(squares), j = 3 (up triangles) j = 4 (down triangles), 
N = 2^2■*, . . . ,2l^ for d = 1 with c = 8 x IQ-^ (hollow 
symbols at the top of the figure), d — 2 with c = 4 x 10"* 
(filled symbols) and d = 3 with c = 4 x 10~^ (symbols with a 
bar at the bottom of the figure). The simulation results are 
averaged over IC configurations for d = 1 and d = 2, and 
over 10* configurations for d = 3. The lines represent the 
(main order) asymptotic theoretical results, namely, \/5 for 
d = 1, 1 for d = 2 and 0.678968 • ■ • for d = 3. 



In kN — ' n 

n=l 



(36) 



and K = I/a/tt 



VII. SUMMARY AND OPEN PROBLEMS 



C. One-dimensional case: rigorous results 

For the one-dimensional lattice one can obtain rigorous 
order-statistics asymptotic expressions for the trapping 
times (lifetimes) from the knowledge of the order statis- 
tics of the diffusion process in the presence of two fixed 
traps | i2L |. without resorting to the Rosenstock approxi- 
mation. We only quote here the results up to first-order 
corrective terms for large N (second-order terms were 
calculated in pl|): 



with Tj^N{m) — Ti_Ar(m) -|- 5j^N{m), 

O 



, ^ ilnlnKA^-7 ^^/In^luKiV 

T\^N{m) = 1 + m- 



In kN 



In kTV 



(34) 



, (35) 



In this review we have discussed some recent advances 
in the field of multiparticle independent random walks. 
The random walk is one of the simplest non-equilibrium 
models in statistical physics, and has played an impor- 
tant role in the study of transport troughout the past 
century Q, |^ |^ . Problems concerning a single random 
walker have been the subject of thorough study, and it is 
today a textbook discipline j^. However, multiparticle 
random walk problems have been posed only very re- 
cently [aQilallfl ^-nd we are just starting to understand 
what happens in this case. We have here dealt exclusively 
with independent random walkers all starting from the 
same site on Euclidean, deterministic fractal and disor- 
dered lattices. Three problems have been discussed: 

(a) The territory problem: the estimate of the average 
number (S-^it)) of distinct sites visited at time t by iV 
random walkers starting from the origin at i = 0. 

(b) Order statistics of exit times from a (hy- 
per)spherical region: the evaluation of the moments of 
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the time tj^N taken by the j-th, j = 1, 2, . . . , random 
walker out of N to arrive at a closed boundary at a given 
distance from the origin. 

(c) Order statistics of trapping times: the estimate of 
the moments of the elapsed time tj^N until the first j 
particles from a set of N are trapped in a disordered 
configuration of trapping sites. 

These problems can not be solved by simple general- 
izations of the solution of the iV = 1 case, and specific 
techniques hava had to be developed for the multipar- 
ticle problems even for independent random walkers. It 
is remarkable that these techniques lead to solutions of 
the aforementioned problems in terms of asymptotic se- 
ries for larg e N which share the same mathematical form 
El III EJEl S 0, i3 . The corrective terms of these 
series decay logarithmically in the number N of random 
walkers, and are consequently, very important even when 
N is very large. This is reflected in the large magnitude 
of the fluctuations on the diffusion front (see Sec. IIV i^ . 

Despite the success in the calculation of the territory 
explored (c.f. Sec. IIV|I . the order statistics of the exit 
times (c.f. Sec.0) and the order statistics of the trapping 
time (c.f. Sec. IVlll . there are still many open problems. 
For example, there exists the problem of the evaluation of 
the distribution of the territory explored or, equivalently, 
the moments {S'^{t)), to = 2, 3, • • • Only the moments 
of the number of distinct sites visited in a given direction 
of a one-dimensional lattice up to time t have been calcu- 



lated rigorously '20|. There are good reasons to believe 
that these moments may be expressed by series in IniV 
of the same form as those in Eg. 0, but this is only a 
conjecture for general media |2l| . The development of 
a procedure for the calculation of these moments in Eu- 
clidean and fractal lattices is a still a challenge. The work 
discussed in this review could also be extended in other 
directions. For example, the extension of our work to 
the case in which the traps are randomly distributed on a 
fractal (deterministic or stochastic) is an interesting open 
problem. Another front opened recently is the study of 
unknotting in granular chains on a vibrating plate [47l |. 
The time corresponding to the unknotting of the first 
knot in a chain with N knots corresponds to the first 
passage time of a random walker of a set of N random 
walkers interacting through an excluded volume pair po- 
tential. This is a good motivation, among many others, 
to investigate the influence of interactions in the quanti- 
ties we have considered in this review, namely, (5Ar(t)), 
(t^^) and (t^^). 
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